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The last decade has seen numerous efforts to achieve imaging resolution beyond that of the 
Abbe-Rayleigh diffraction limit. The main direction of research aiming to break this limit seeks to 
exploit the evanescent components containing fine detail of the electromagnetic field distribution 
at the immediate proximity of the object. Here we propose a solution that removes the need for 
evanescent fields. The object being imaged or stimulated with sub- wavelength accuracy does not 
need to be in the immediate proximity of the super lens or field concentrator: an optical mask can 
be designed that creates constructive interference of waves known as superoscillation, leading to a 
sub- wavelength focus of prescribed size and shape in a 'field of view' beyond the evanescent fields, 
when illuminated by a monochromatic wave. Moreover, we demonstrate that such a mask may be 
used not only as a focusing device, but also as a super-resolution imaging device. 



The last decade has seen numerous efforts to achieve 
imaging resolution beyond that of the Abbe-Rayleigh 
diffraction limit, which proscribes the visualization of fea- 
tures smaller than about half of the wavelength of light 
with optical instruments. The main direction of research 
aiming to break this limit seeks to exploit the evanescent 
components containing fine detail of the electromagnetic 
field distribution. Indeed, many powerful concepts like 




scanning near-field optical microscopy (SNOM) 111, 
use of various forms of field concentrators 0, 0, [_ 
superlenses @, H, 0, QU EH, EH, EH depend on the ex- 
ploitation or recovery of evanescent fields. The proper 
far-field optical superlens Q requires bulk negative in- 
dex materials that are still to be developed, while other 
designs, though offering substantial advances, are united 
by a common severe limitation: that the object being 
imaged or stimulated must be in the immediate proxim- 
ity of the superlens or field concentrator. Here we pro- 
pose a solution that removes this limitation: an optical 
mask can be designed that creates constructive interfer- 
ence of waves known as superoscillation, leading to a sub- 
wavelength focus of prescribed size and shape in a 'field 
of view' beyond the evanescent fields, when illuminated 
by a monochromatic wave. Moreover, we demonstrate 
that such a mask may be used not only as a focusing 
device, but also as a super-resolution imaging device. 

In fact, in his seminal 1952 paper Toraldo Di Fran- 
cia demonstrated that propagating waves can create sub- 
wavelength localization of light in the far-field with the 
suggestion of a pupil design providing an accurately tai- 
lored sub-wavelength diffraction spot using a series of 
concentric apertures [14]. Such a sub- wavelength con- 
centrator could be employed as a focusing device in a 
super-resolution scanning optical microscope where the 
object is placed several wavelengths away from the de- 
vice, thus removing the main limitation of near-field in- 
struments. More recently, Berry and Popescu [l5[, start- 
ing from earlier works on quantum mechanics, showed 
that diffraction on a grating structure could create sub- 
wavelength localizations of light that propagate further 
into the far field than more familiar evanescent waves. 



They relate this effect to the fact that band-limited func- 
tions are able to oscillate arbitrarily faster than the high- 
est Fourier components they contain, a phenomenon now 
known as superoscillation [l6|, E3, E3 . Examples of sub- 
wavelength localizations of light generated by a nano-hole 
array pj| [20[ and a thin meta-dielectric shell [2l| have 
been demonstrated recently. Research on beating the 
diffraction limit actually has an even longer history: in 
1922, Oseen, with reference to Einstein's radiation 'nee- 
dle stick', proved that a substantial fraction of the emit- 
ted electromagnetic energy can be sent into an arbitrar- 
ily small solid angle 1221 ] . Beginning from the pioneering 
work of Shelkunoff [23], the microwave community con- 
templated the idea of achieving antennae that beat the 
diffraction limit for directivity: several authors were able 
to prove that for a linear array of properly adjusted ra- 
diating antenna dipoles, there were no theoretical limits 
to directivity whatsoever [24|, [25[. However, the sharp 
increase in the proportion of reactive to radiated power 
that would be required to achieve super-directivity means 
that the antenna gain increase is offset by the need to pro- 
vide an even higher increase in the power to the antenna 
to maintain the signal level, thus rendering the concept 
of super-directive antennae impractical. As we will see 
below, achieving a sub- wavelength localization of light in 
the far-field also comes at a price of losing most of the 
optical energy into diffuse sidebands. Nevertheless, op- 
tical microscopy applications can tolerate much higher 
losses than those acceptable in antenna design: scanning 
microscopes can work with only a few photons per sec- 
ond, giving one around 19 orders of magnitude of power 
reserve (assuming that a 1 W laser is used as the optical 
source). 

In this letter we make a further step in the investiga- 
tion of the potential of superoscillation for imaging and 
achieving sub-wavelength foci. We derive an algorithm 
for designing a mask that creates a sub- wavelength focus 
of prescribed size and shape within a prescribed 'field of 
view' when illuminated by a monochromatic wave. More- 
over, we show that such a mask may be used not only as 
a focusing device, but also as a super-resolution imaging 
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FIG. 1: Optical superoscillation. An arbitrary field distribu- 
tion with an arbitrarily small hot-spot 5x, within a limited 
area [—D/2, D/2], can be generated by diffraction of a plane 
electromagnetic wave from a purposely designed mask. 



tool. We also study the role of manufacturing imper- 
fections on the achievable super-resolution and suggest a 
design for a superoscillation plasmonic energy concentra- 
tor. 

The typical situation that we address here is presented 
in FigHJ we aim to design a mask which, within a lim- 
ited area [—D/2, D/2] (field of view), will create a small 
hot-spot of light concentration (superoscillation) with a 
width Sx located outside the evanescent zone, at a dis- 
tance z > A from the mask, where A is the wavelength of 
light illuminating the mask. We argue that in principle 
a mask can be designed to create a hot-spot that is ar- 
bitrarily small, with an arbitrary profile, located at any 
given distance from the mask. 

Our consideration is limited to a one-dimensional mask 
t(x) creating a one-dimensional sub- wavelength field dis- 
tribution h(x) (superoscillation) when illuminated by 
plane wave at wavelength A, although generalization to 
a two-dimensional case is trivial. The desired superoscil- 
lation feature located at a distance from the near-field 
zone of the mask can only be created by diffraction on 
the mask if the feature can be decomposed into a se- 
ries of plane waves with wave- vector |kn| = There- 
fore, the main step in designing the mask is to present 
the desired superoscillating field as a series of bandwidth 
limited functions that can be decomposed into free-space 
plane waves of the given wavelength A. We argue that 
any arbitrary small field feature can be presented as a 
series of band-limited functions if we are concerned with 
a prescribed 'field of view' [—D/2, D/2]. This may be 
achieved using the formulism of prolate spheroidal wave 
functions developed by Slepian and Pollark [26[ to treat 
problems of information compression. This is a complete 
set of functions orthogonal in the interval [—D/2, D/2] 
and across the whole range [—00,00]. The main feature 
of prolate spheroidal wave functions is that they are band- 



FIG. 2: Normalized intensity profiles for a hot-spot with Sx = 
0.07A (black dotted line) and its approximations by series of 
bandwidth- limited prolate spheroidal wave functions \h(x)\ 2 
truncated at: N — 2 (red line, £x2=0.38A); N = 6 (blue line, 
£x6=0. 21A); N = 10 (pink line, £xi =0.14A); N = 26 (green 
line, 5x26 — 0.07A ). The dashed line shows the A/2 limit for a 
hot-spot achievable with a high-numerical aperture cylindrical 
lens. Insert: Energy (£) contained within the hot-spot as a 
function of its size (Sx) 

limited to a frequency domain [— fco?fco]- Therefore the 
mask design algorithm comprises the following steps: ini- 
tially the desired sub-wavelength hot-spot is presented 
as a series of prolate spheroidal wave functions, which 
can be truncated when a satisfactory level of approxima- 
tion is achieved; at the second step this series of prolate 
spheroidal wave functions is presented as a series of plane 
waves and, using the scalar angular spectrum descrip- 
tion of light propagating from the mask to the super- 
oscillating feature, the required complex mask transmis- 
sion function t(x) can be readily derived. The formulism 
of the algorithm is presented in the Appendix. 

In what follows, we will give an explicit example of 
a mask designed to generate a sub-wavelength concen- 
tration of light. Let us aim for a single hot-spot field 
distribution Sinc(ax/X) centered in the 'field of view' 
[-D/2,D/2]. The full width at half maximum Sx of 
the intensity profile of this distribution is measured as 
2. 784 A/a. From now on we will set a = 40, aiming there- 
fore to achieve a Sx = 0.07A hot spot in a D = 1.2 A 
field of view, far beyond the Abbe- Ray leigh limit. Fig. 2 
shows the intensity profile of this distribution alongside 
a number of consecutive approximations to the distribu- 
tion formed by limited series of prolate spheroidal wave 
functions, and a curve representing the Abbe- Ray leigh 
limit that would be achievable by a high-numerical aper- 
ture cylindrical lens. One can see that the series rapidly 
converges and that when N = 26 the width of the ap- 
proximation is practically the same as that of the target 
field distribution. 

Fig. 3 shows the intensity |t(x)| 2 and phase cp profiles 
of the mask t(x) = \t(x)\e t(p required to create the field 
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FIG. 3: Mask profile for the generation of a sub- wavelength 
hot-spot, (a) shows intensity |t(x)| 2 (blue line) and phase 
<p(x) (red lines) profiles of the mask transmission function 
t(x), which generates a hot-spot with Sx = 0.21 A at a distance 
z = 20A from the grating. At the bottom, a grey-scale map 
of the intensity profile is also presented, (b) shows close-up 
detail of kinks in the phase curve at the points highlighted 
with arrows in (a). 
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profile corresponding to N = 6, which has Sx = 0.21A, at 
a distance z = 20A from the grating. As the mask trans- 
mission function is even, only part of profile for x > 
is shown here. One can see that there is a small cen- 
tral area [— 40A, 40A] of the mask that transmits most of 
light. Beyond that area, there is a low-intensity broad 
shoulder of a slowly fading transmission characteristic. 
The phase profile of the mask resembles that of a con- 
cave lens where the optical thickness is at a minimum in 
the centre. The monotonous increase is only interrupted 
at a few positions where transmission amplitude is close 
to zero (indicated by arrows) and the phase shows a kink 
as illustrated in the zoomed sections (b). 

The superoscillation process is based on the precise 
and delicate interference of waves, so we investigated its 
stability with respect to the manufacturing tolerances to 
which masks can be fabricated and to the need to use 
masks of finite length (see the description of the practi- 
cal device design below). We found that for the mask 
presented in Fig. 3 an 8000A long device may be used in- 
stead of an infinitely long mask without any substantial 
degradation of performance (Sx= 0.21A for both infinite 
and truncated masks with a 5% intensity decrease for the 
finite mask). A practical mask is likely to be manufac- 
tured by electron-beam lithography with a certain lim- 
ited resolution, so we also investigated the dependence 
of superoscillation hot-spot size and intensity on the pix- 
elation of the grating design: we replaced the smooth 



FIG. 4: The role of finite mask length and manufacturing 
imperfections on mask performance, (a) shows the intensity 
and phase characteristics of an ideal mask (smooth line) and 
its stepped equivalent, where the step of 0.16 A mimics the 
manufacturing pixelation. (b) Calculated intensity profiles of 
hot-spots generated by an ideal infinitely long mask (black 
line); an ideal mask truncated to a length of 8000A (red line); 
and 8000 A long masks pixelated to 0.08 A (cyan line), 0.1 2 A 
(blue line), 0.16A (green line ) and 0.2 A (pink line). 



transmission function with a stepped functions in which 
the step width corresponded to the resolution P of the 
manufacturing process (see Fig. 4a). Fig. 4b shows the 
variation of the intensity profile with increasing pixela- 
tion. The superoscillation process is remarkably stable 
against this manufacturing imperfection: with pixel size 
increasing to 0.2 A, the Sx width of the hot-spot only in- 
creases by about 8%. Pixelation has a more substantial 
influence on the peak intensity, which is reduced by about 
40% under such conditions. 

To quantitatively characterize the optical energy con- 
tained in a superoscillating hot-spot, we define £ as the 
ratio between the energy contained in the hot-spot (be- 
tween the first intensity minima either side of the cen- 
tral peak) and the total energy transmitted through the 
mask. In the inset to Figj2] we show values of £ against 
the width of the peak Sx for different levels of approxi- 
mation N = 2, 6, 10. One can clearly see that the smaller 
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FIG. 5: (a) Superoscillation mask as a sub- wavelength imag- 
ing device, (b) Intensity profile of the image of two sub- 
wavelength slits separated by a distance Ax± = 0.24 A, located 
at a distance 20 A from the mask. The intensity distribution 
is also shown in gray-scale at the bottom of the image where 
the dotted lines indicate the positions of the two slits. For 
comparison, the dashed line shows the intensity profile of the 
image field of the two slits as seen imaged by a conventional 
cylindrical lens with unitary numerical aperture and a reso- 
lution of A/2. 



the hot-spot 'pedestal' (the background level within the 
field of view), for instance by stipulating only the maxi- 
mum field intensity in the pedestal area relative to that 
of the hot-spot, but not the exact profile of the pedestal. 

Recently it was shown that a hole array can be used 
as an imaging device [27]. Similarly, a super-oscillating 
mask can also be designed to image a sub-wavelength 
object. This is illustrated in Fig. 5 where we show a 
mask with a transmission function that converts an ob- 
ject located 20 A from the mask to an image on the op- 
posite side of the mask, also at 20A (see formula A6 of 
the Appendix). Fig. 5(b) illustrates the imaging of two 
incoherent slit sources (each 0.04A wide) separated by a 
distance of 0.24A. Here the slits are clearly resolved ac- 
cording to the Rayleigh criterion [28], which state that 
the total intensity at the saddle point of the sum inten- 
sity profile of two just-resolved slit sources is 81% of the 
maximum intensity. 

Regarding practical implementations of the superoscil- 
lating mask, the above calculations demonstrate that 
manufacturing tolerance should be of the order of 0.1 A. 
Thus, manufacturing such a phase mask from a slab of di- 
electric material for microwave and THz frequency focus- 
ing should not be a very challenging problem. It can then 
be covered with an absorbing film of variable thickness 
to create the desired transmission profile. Fabricating 
such a mask for the IR and optical parts of the spectrum 
presents a more significant challenge and a fabrication 
accuracy of between 5 and 50nm will be required. How- 
ever, even such challenging phase masks may be created 
from glass with the diamond milling techniques used to 
produce aspheric lenses, while further fine-tuning with 
true nanoscale resolution may be achieved using focusing 
ion beam milling. Here the absorbing part of the mask 
could be a metal film of variable density prepared by UV 
or e-beam lithography. 




the hot-spot is, the lower the proportion of energy that 
goes into it. For example, for N = 6 only 1.74 x 10 ~ 5 of 
the total transmitted energy will be focused into the hot- 
spot. The energy contained in the hot-spot is also signifi- 
cantly dependent on the field of view: if the required field 
of view is decreased by a factor of two from D = 1.2 A 
to D = 0.6A, the proportion of energy going into the 
hot-spot at TV = 6 increases by nearly three orders of 
magnitude to 4.8 x 10 -2 . In fact, a more general consid- 
eration of the super-oscillating functions shows that the 
intensity in the hot-spot may only decrease polynomially 
with its width [lTj • We therefore argue that when design- 
ing practical realizations of super-oscillating masks, less 
demanding requirements on the shape of hot-spots may 
yield much higher values of £. In particular, one may 
achieve much higher values of £ for a given size of hot- 
spot by relaxing the requirements on the field structure of 




FIG. 6: Plasmonic superoscillation focusing. A plasmonic 
mask is archived by placing a sculptured dielectric refractive 
layer on top of a flat metal surface to create the phase profile 
followed by a sculptured lossy layer to create the intensity 
profile. 

The idea of superoscillation could also be applied to 
the creation of a super-high-resolution plasmonic device. 
Surface plasmon-polaritons are collective oscillations of 
light and electrons that propagate along the interface be- 
tween a metal (e.g. gold, silver) and a dielectric. They 
are essentially two-dimensional waves that can create 
complex field patterns by interference [29|, [3(j • As plas- 
mon wavelengths are somewhat shorter than those of 
light at the same frequency, plasmonic devices promise 
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better resolution than optical devices even within the 
conventional Abbe- Ray leigh diffraction limit. We argue 
that this may be enhanced even further through the use 
of superoscillation as the approach developed above may 
be easily applied to plasmons. A possible implementa- 
tion of a plasmonic focusing device is presented in Fig. 6. 
Here, in order to create a desired intensity and phase 
profile of the mask we exploit the fact that the complex 
refractive index for plasmons depends on the dielectric 
properties e of the medium forming an interface with the 

metal (3l| : n = y/ / m _^ £ . By preparing a film of lossless 

dielectric on a metal surface with a prescribed profile, 
one can create the desired phase profile for the plasmon 
mask. Using an additional lossy profiled dielectric layer 
on the top of the first the necessary intensity profile can 
be engineered. Here the transmission function of the sur- 



face plasmon mask should be derived taking account of 
losses (see Appendix, formula A7). 

In summary, we have shown that an optical mask can 
be designed that creates a a sub-wavelength focus in an 
area beyond the evanescent fields. Such a mask may 
also be used as a super-resolution imaging device with 
numerous applications for instance for imaging inside a 
leaving cell which is impossible with a near-field device. 
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where i?g^(c, 1) is a radial prolate spheroidal wave 
function of the first kind and 



/«0 
H N (u)e iux du (A4) 

Comparing formulae (|Al| and (|A4j) , we find that the 
required transmission function t(x) of the mask is 



APPENDIX A: DESIGNING A 
SUPEROSCILLATING MASK 

Here we describe a method to design a mask with a 
complex transmission function t{x) that will generate a 
prescribed field distribution f(x) within a limited region 
[—D/2,D/2] at a distance z from the mask using pro- 
late spheroidal wave functions [26]. We assume that the 
mask is illuminated at normal incidence with a plane 
monochromatic wave E(x,z = 0) = 1 (the time depen- 
dent factor e tUJt is omitted) with a wavelength A = 2tt //co- 
in the scalar angular spectrum description of light prop- 
agation [32[ the field at a point (x, z) is 



E(x,z)= f° T(u) 

J — k 



du 



(Al) 



where T(u) is the Fourier transform of t(x). We now 
approximate h(x) = E(x,z) as a limited series of orthog- 
onal prolate spheroidal wave functions %jj n (c,x) that are 
band-limited to the frequency domain [— fco? &o] 



n—N 



h N {x) = ^ CL n (c)^n(c, X) 
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Here a n depends on a constant c = 2jP and the Fourier 
transform function of hjy(x) is given by 



n—N rk 
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It follows from the angular spectrum representation 
that to design a mask with a transmission function m(x), 
which upon illumination by a single slit source (located at 
a distance z\ from the mask), will convert its divergent 
incident field E (x) into a prescribed field distribution 
h(x) at a distance z from the other side of the mask, the 
following formula should be used: 



m(x) 



t(x) 



f*° ko F(u)e iux e izl V k2 o- u2 du 



(A6) 



where F(u) is the Fourier transform function of the slit 
source. 

In the case of lossy medium the transmission function 
A(5) shall be corrected: 



t(x) 



r\ko\ iranil> n (c,^ 

x e iux-izRe[^ kl~u 2 ] ^ u 



(A7) 



where fco is a complex wave- vector in the lossy medium. 



